Spin-charge interplay in electronic liquid crystals: fluctuating spin stripe driven by 

charge nematic 
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We study the interplay between charge and spin ordering in electronic liquid crystalline states with 
a particular emphasis on fluctuating spin stripe phenomena observed in recent neutron scattering 
experiments [J, 0] ■ Based on a phenomenological model, we propose that charge nematic ordering 
is indeed behind the formation of temperature dependent incommensurate inelastic peaks near 
wavevector (7r,7r) in the dynamic structure factor of YBa2Cu306+ H . We strengthen this claim by 
providing a compelling fit to the experimental data which cannot be reproduced by a number of 
other ordering possibilities. 

PACS numbers: 71.10.Hf,71.45.Gm,75.10.Jm 
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Kaleidoscopic variety of competing ordering tendencies 
is both a hallmark of correlated electron fluids, such as 
cuprate and Fe based superconductors d, H|, and a theo- 
retical challenge. As such, a clear identification of a bro- 
ken symmetry phase offers a valuable guiding principle. 
Recent observations of temperature, energy and doping 
dependent onset of anisotropy in inelastic neutron scat- 
tering (INS) studied by V. Hinkov et al. [l[ provides an 
opportunity for just such identification. 

The symmetry of the "fluctuating spin stripe" phenom- 
ena (one-dimensional incommensurate spin modulation 
at finite energy) observed in Ref. [l| is consistent with 
that of a nematic phase [B| (a metallic state that breaks 
rotational symmetry without breaking translational sym- 
metry). Furthermore, the qualitative departure in the 
magnetic response of underdoped YBa2Cu30e+yin Ref. 
l[ from that of optimally and overdoped regimes y > 0.5 
6H2| indicates the possible existence of a quantum crit- 
ical point at around y ~ 0.5 as we sketch in Fig. [T] 
In specific, low-energy features are enhanced in INS of 
y = 0.45 while the extensively studied high energy "hour- 
glass" dispersion which are prominent at higher doping 
p42( is suppressed. (Considerable attention has been di- 
rected towards this resonance feature and its significance. 
See e.g. Ref.fiol]). However, despite the reported temper- 
ature dependence of the finite-frequency incommensura- 
bility being suggestive of an order parameter [111 ], it has 
not been clear how this quantity can be related to a spe- 
cific order parameter since an order parameter is defined 
by broken symmetry of the ground state. 

There has been a number of theoretical studies re- 
garding possible signatures of electronic liquid crystal 
physics in the magnetic response 12Ul4j]. While these 
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studies shed light on the hour-glass dispersion observed 
in y > 0.5 at high energies, their connection with the 
low-energy phenomena in the underdoped regime with 



FIG. 1: (color online) Schematic phase diagram of 
YBa2Cu306+ H . Here the shaded region bound by Tjv repre- 
sents a nematic phase and the dashed line below Tn represents 
the crossover temperature T* to a "fluctuating spin stripe" 
behavior. Sufficiently strong nematic ordering at low enough 
temperatures could further stabilize a static spin stripe phase. 



y < 0.5 is unclear. Moreover, they focused on the super- 
conducting phase while the observed onset of fluctuating 
spin stripe behavior is at Tn ~ 150K, well above the 
superconducting ordering temperature T c — 35K. 

In this letter we propose that charge nematic ordering 
is the driving force behind the fluctuating spin stripe phe- 
nomena observed in underdoped YBa2Cu3C>6+y . We con- 
sider a metallic system proximate to antiferromagnetic 
(AFM) ordering and show that charge nematic ordering 
quite uniquely can induce fluctuating and even static spin 
stripes thus providing a concrete connection between the 
charge [l5j and spin aspects of liquid crystalline behavior 
in underdoped YBa2Cu3C"6+ y . Our claims are supple- 
mented by a successful fit with available INS data[l|. 

Phenomenological Model We start by noting that V. 
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Hinkov et al. [l| detect the dynamic onset of anisotropy 
through incommensurate inelastic peaks near the AFM 
wave vector Q = (tt, 7r) in a metallic system (see FigJ5]). 
Given a microscopic theory of itinerant magnetism being 
an open question, we take a phenomenological approach 
as a first pass through the problem. In the presence of 
long-range AFM ordering, low-energy excitation in the 
particle-hole channel is dominated by gapless spin waves 
near Q with infinite life time: 
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^(k+Q+q,0+a;)o ! Q/J V/3(k,fi). (1) 



Here the operators V'KV'a) are the fermion creation (an- 
nihilation) operators with spin a =f, J,, tfi and a are 
vectors in spin space with a l aj3 representing (a, (3) com- 
ponent of Pauli matrix a % for i = x, y, z. q denotes the 
wavevector of the spinwaves measured from Q and uj de- 
notes their frequency. For underdoped cuprates outside 
of but close to the AFM phase (see FigH]), this spin- 
wave can be either damped or gapped and well defined. 
To quadratic order in 0, the appropriate effective action 
takes the form 16 [ 
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where g is an overall scale and A (q) is the momentum 
q dependent spin-wave gap in the absence of damping 
and r is the damping energy scale. Here we take the 
damping to be of Landau damping form 17[ but most of 
our conclusions below are insensitive to the details of the 
damping dynamics. 

The dynamic structure factor measured by INS exper- 
iments will be proportional to the spectral function 
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for the action of Eq. @. The uniform component of the 
gap A(q = 0) sets the energy scale above which INS in- 
tensity is significant and the q dependence of A(q) deter- 
mine the distribution of intensity in the Fourier space. In 
the absence of other symmetry breaking tendencies, A(q) 
will be a minimum at q = and respect the point group 
symmetry of the system. Hence the INS intensity above 
A(q = 0) will be C4 symmetric (with static orthrhom- 
bicity due to chain layer) and peaked at q = 0. Tem- 
perature and energy dependent anisotropic incommensu- 
rability observed in Ref.|l( indicates additional ordering 
tendencies at play. We first study the effect of charge 
nematic ordering motivated by transport anisotropy ob- 
served in the nearby regime [15|. 

In a nematic fluid, charge degrees of freedom collec- 
tively break rotational symmetry of space while preserv- 
ing the translational symmetry (hence remains metallic) 



and any other symmetry of the system such as time re- 
versal and spin rotation^. A nematic order parame- 
ter in a continuum system has the symmetry of charge 
quadrupole moment (1 = 2 representation of 5*0(2) ro- 
tational group). It has two real components and it can 
either be written as a symmetric traceless tensor of rank 
two or as a complex field [lil. Hij. 

As the nematic ordering involves spatial symmetry 
breaking, the nature of the order parameter itself is af- 
fected by crystal fields due to the lattice which lowers 
the spatial rotational symmetry to a discrete group. For 
instance in a square lattice with C^ v symmetry, the or- 
dering can occur in one of two channels d x 2_ y 2 or d xy 
reducing the order parameter symmetry to Ising-like (a 
single component real field) in either case [ifj. Due to this 
discretization of the order parameter symmetry, the ne- 
matic fluctuation below the transition temperature T/v 
will be massive. In YP^CuaOg+y, the weak external 
field imposed by the chain layer will likely pick the ne- 
matic to occur in the d x 2_ y 2 channel and a representative 
form of order parameter can therefore be written as 



N = 
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N(cosk x - cosfc g )'0j,(k)'0a(k). (4) 



However, it is important to note that any electronic quan- 
tity that is odd under 90° spatial rotation such as ef- 
fective mass anisotropy ratio (m y — m x )/(m y + m x ) or 
transport anisotropy (p xx — Pyy)/(pxx + Pyv} can serve 



as the order parameter for the nematic state 18 1. In fact, 
short of probes that couple directly to charge quadrupole 
moments, nematic phases so far has been mostly de- 
tected through temperature dependent in-plane trans- 
port anisotropy in quasi 2D systems 2l|, \22 \ . 

In YBa2Cu306+y, Ando et al. [HI observed that trans- 
port anisotropy increases upon under doping below y ~ 
0.5 while the effect of CuO chains diminishes. Its elec- 
tronic origin is further supported by the fact that this 
anisotropy only sets in below a doping dependent onset 
temperatures TV (180K for y = 0.45) making nematic 
ordering a compelling candidate for background broken 
symmetry in underdoped YBa2Cu30g+ !/ as advanced in 
Ref. [l[ without a theoretical model. 

Consider the effect of d x 2_ y 2 nematic ordering below 
Tjv on the spin fluctuations in Fig. [T] The effective action 
now depends on both (f> and N . Away from the classical 
critical region near T/v, the gapped nematic fluctuations 
can be integrated out and the effect of finite TV can be 
represented by ./V-dependent coefficients in the gradient 
expansion of S[cf>, N] in the long distance limit. We rep- 
resent such N dependence of S[<j>, N] using the N depen- 
dent "gap" A(q;iV). On symmetry grounds A(q;7V) in 
the long distance limit takes the following form to the 
quadratic order in q: 



A 2 (q;N)=Al(N)+c 2 (N) q 2 ~c 2 2 (N)N( q l- q 2 y ) 



(5) 
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where all momenta are in units of the lattice constant (i.e. 
q x = k x a). Assuming tetragonal symmetry of underlying 
lattice in the absence of symmetry breaking field, all the 
functions A , Cq , c 2 should be even functions of A since 
N — > —A under 90° spatial rotation. For small A, we 
expand A(q; A) in powers of N and treat A ,co,c 2 as 
independent of N to lowest order in TV [23| . (Note that 
this approach would be only valid up to the nematic gap 
scale which we estimate to be about ~ 18 meV from the 
onset temperature for transport anisotropy TV ~ 180 K.) 

Notice that the nematic ordering allows for anisotropic 
flattening of the momentum dependence of A(q; N) in 
Eq. © which will elongate the inelastic (tt, it) peak in the 
x direction. For N > N* , with A* = c 2 /c|, this effect 
will further shift the peaks to incommensurate positions 
at (n ± S,tt) at low-energy (w < A(q) for any q) due to 
a dip in A 2 (q; N) at q = (±S(T), 0) with 

8(T) oc (N(T) -A*) 1/2 . (6) 

Such incommensurability will naturally have tempera- 
ture dependence resulting from the temperature depen- 
dence of N(T). We represent this crossover temper- 
ature T* for onset of such fluctuating spin stripe at 
N(T*) — A* as a dashed line in Fig. [U It is noteworthy 
that while the form of Eq. ([6]) is reminiscent of an order 
parameter, there is no singularity in the free energy at 
T* (this is a crossover). 

At even lower temperatures, large enough A may 
even stabilize static spin stripes. While our lowest or- 
der in A analysis may not be applicable in its explicit 
form for large A, it predicts a critical value N c where 
A(±S, 0; N c ) = 0. For A > A c the system will develop 
static spin stripes. Such static spin stripes were recently 
observed at T = 2 K by D. Haug et al. Q and a similar 
effect was also reported in La2- x Sr x Cu04|24] 
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FIG. 2: (color online) Fit to momentum dependence of INS 
intensity along a (blue dots) and b (red triangles) axes. Data 
obtained from Ref.[lJ taken at T = 5K and lu — 3meV. Note: 
any slight asymmetry observed in the data is not accounted 
for in the phenomenological model presented here. 

Comparison with experiments The model defined by 
Eqs. © and © provides a natural connection between 
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FIG. 3: (color online) The frequency dependence of the mo- 
mentum integrated spectral function. (Note: this is really only 
valid up to the nematic fluctuation mass scale of 18 meV). 

key observations of transport measurements [l5| and re- 
cent INS studies in YBa2Cu306+y One non-trivial 
prediction is the temperature dependence of the incom- 
mensurability 5(T) in the vicinity of T* . 

It is straightforward to show that Eq. © implies tem- 
perature dependence of 5(T) reminicsent of that of a 
mean field order parameter yet without any thermody- 
namic singularity at T — T* . This is independent of the 
unknown order parameter exponent j3 of N oc (T — T^)^ 
and any microscopic details. As long as N* is small 
and the spin waves are gapped, so that we may use the 
quadratic action of Eq.((2]), the (N — A*) 1 / 2 dependence 
of S in Eq.© is valid to leading order in (A — A*). By 
inserting N oc (T — T^) 13 to Eq.© , one can show 

S oc (T-T*) 1 / 2 . (7) 

which is in good agreement with experimental observa- 
tion in Ref.[l|. Thus we have established a model which 
yields observed T-dependent incommensurability. 

As a practical test of ourphenomenological theory, 
we fit the INS data of Ref.fl] with the spectral func- 
tion Eq.© using A(q; A) from Eq.©. Figj2] shows a fit 
to momentum line cuts of INS data at u — 3 meV and 
T = 5K by setting g/T = 5.1 meV, [A 2 - (3meV") 2 ]/r = 
2.7 meV, cl/T = 16 meV, cjA/r = 24 meV. (We 
also include symmetry allowed quartic terms of the form 
M«2 + ^qlql + A 3 A(g4 - q *) to fit high q part 

of data [25|.) Note that this momentum dependence fit 
at oj — 3meV is insensitive to the relative strength be- 
tween damping and gap energy scales T/Aq. Only the 
frequency dependence is sensitive to this ratio. Similar 
fits at higher temperatures (T = 40K and T = 100K) hint 
at thermal fluctuation driven damping at higher tempera- 
tures and T* being below 100K (see supplementary online 
material for details). Note that we are using a different 
scheme for extracting T* from that of Ref. [l( . 

Fixing most of model parameters by the u = 3meV, 
T = 5K data, we fit the frequency dependence of uni- 
form susceptibility x"( w ) a t 5K. The energy dependence 
fit shown in FigJ3] is valid up to the estimated scale 
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of nematic fluctuation mass (« 18 meV). Both FigsfS] 
and [3] show good agreement with data. However given 
quite a few parameters, we should stress that fitting is 
more of a check rather than the main result of this pa- 
per. Nonetheless, the frequency dependence allows us 
to estimate T/Ao ~ 1.3 giving a reasonable fit. Note 
that while one expects the spinwave to be well defined 
inside superconductor due to lack of low-energy excita- 
tions to scatter off, the data of Ref.[l[ appears to in- 
dicate significant amount of damping. Speculating the 
source of this damping is beyond the scope of our pa- 
per but this might be consistent with existence of the 
low-energy excitations reported in superconducting un- 
derdoped Bi 2 Sr 2 CaCu2O 8+1 5 0. 

Effect of Other Ordering Tendencies For complete- 
ness, we also consider the effect of two other competing 
orders possibly proximate to the regime of our interest: 
the charge stripe(smectic) and d-wave superconductivity. 

For charge stripe ordering with wave vector Q cs the 
order parameter is 

pq. s = J ^{^(k + Q M )^(k)), (8) 

where the spin index a is summed over. We therefore 
find an additional contribution to the spinwave action 




(9) 

From Eq. (J9)) it is clear that the existence of charge stripe 
shifts the spinwave gap by a constant (g cs term) and it 
can also induce fluctuating and static spin stripe with in- 
commensurability fixed by the charge stripe wave vector 
by S = Qcs/2 {g'cs term). However, incommensurability 
so induced will be T independent since it is fixed by Q cs . 

In the presence of d-wave pairing the leading addi- 
tional contribution to the spinwave action is 

SS sc = -g sc \A d f J ^|0( q ,u;)| 2 . (10) 

Hence while d-wave superconductivity can also shifts the 
commensurate spinwave gap A(q = 0) by a constant 
g sc \Ad\ 2 and enhance (for g sc < 0) or suppress (for 
g sc > 0) it cannot lead to the observed T dependent in- 
commensurability either. These conclusions can be gen- 
eralized to other pairing symmetries. 

Discussions We have shown that the existence of 
charge nematic ordering can take a system proximate 
to AFM and induce fluctuating spin stripe phenom- 
ena: incommensurate INS peaks. Within a symmetry 
based phenomenological approach we explained the ob- 
served (T* — T) 1 / 2 dependence of the incommensurabil- 
ity. Hence, we established a concrete model connecting 



the spin and charge aspects of the electronic liquid crys- 
talline behavior of YBa2Cu306+ y . We further argued 
that this is a unique feature of charge nematic ordering. 

One natural implication of our model is the possibil- 
ity that outside but close to the charge nematic phase, 
impurities could stabilize fluctuating or even static spin 
stripes. Since impurities break the crystal lattice sym- 
metry, they naturally provide a symmetry breaking field 
which will induce a finite N. In weak coupling, this would 
be given by N — xn^n where /ijy captures the orienta- 
tional symmetry breaking of the impurities (see Ref. 5 
for a similar argument applied to charge stripes). Since 
the response xn to this perturbation is expected to be 
large near a nematic phase, an N larger than N* could 
be induced and fluctuating spin stripes will result. If it is 
further larger than N c , static spin stripes will be induced. 

Several future directions include, establishing a mi- 
croscopic model which reproduces our phenomenological 
model in the long-wavelength limit and studying the ef- 
fects of a magnetic field. This latter question is particu- 
larly important in the context of observations of a mag- 
netic field stabilizing static spin stripe phase 0. Since 
such a phase is accessible from our model, it will be in- 
teresting to include the effect of magnetic field. 

Furthermore, our analysis may provide a starting point 
for investigating the interplay between spin and elec- 
tronic liquid crystalline ordering in other systems. By 
now there are a growing number of candidate corre- 
lated systems for such interplay, including Mn-doped 
Sr 3 Ru207[i3] and Fe-based superconductors fijj] . 
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Supplementary Material 

Here we discuss higher temperature fits to the inelastic 
neutron scattering data of Ref.[l[. Short of explicit tem- 
perature dependence in our phenomenological model, we 
consider two approaches: (a) we assume implicit temper- 
ature dependence of all symmetry allowed phenomeno- 
logical parameters used in fitting the low temperature 
(T = 5K) data; (b) we assume the dominant effect of 
temperature is to introduce static damping due to ther- 
mal fluctuation. 

Assuming implicit T- dependence of fitting parame- 
ters, we fit the momentum dependent neutron scattering 
data taken at temperatures T = 40K and 100K. For the 
fits shown in Fig. 21 the parameters are g/T — 5.8 meV 
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(8.3 meV), [Ag - (3meT/) 2 ]/T = 2.9 meV (3.9 meV), 
cl/T = 10 meV (15 meV), cjN/T = 15 meV (5 meV), 
Ai/r = 26 meV (14 meV), and A 3 /r = 18 meV (2 meV) 
for T = 40K ( and for T = 100K respectively). 

At higher temperatures we observe from this fit that 
c > c 2 at T = 100K while c < c 2 for both 40K and 5K. 
As Co < C2 is necessary for a double peak structure sepa- 
rated by S (our definition of incommensurability, see Eq. 
(3) of text), we interpret the 100K fit to be above the 
crossover temperature scale T* below which the struc- 
ture factor shows a double peak. Assuming T)v ~ 150-K", 
the system is then in a nematic phase with anisotropic 
structure factor at T = 100K and above the crossover to 
the low temperature fluctuating stripe behaviour. 



o 



















Q+q 


T 71 T 

Q+q 



FIG. 4: Fit to momentum dependence of INS intensity along 
a (blue dots) and b (red triangles) axes. Data obtained from 
Ref.[l| taken at (a) T = 40K and (b) T = 100K with lo = 3 
meV. 

Consider the second approach to fitting mentioned 
above. Since T = 40K > hw/ks for the neutron energy 
transfer of w = 3meV, sizable damping due to thermal 
fluctuation is expected at 40K compared to T = 5K. 
Hence, we incorporate thermal fluctuation driven damp- 
ing into the effective action through Tq(T) such that 
r (T = 5K) = 0: 

S[ ^ ] = £ / ( ir ° + + ^ " A ' (q) ) l<?(q ' W)|2 - 

It turns out both 5K data and 40K data can be fit by fix- 
ing all other parameters except overall intensity to their 
5K values and setting T (T = 40K)/T = 1.4 meV (see 
Fig. [5]). This signifies enhanced damping as the main 
difference between these data sets . 

Both the experimental resolution and the phenomeno- 
logical nature of our theory make it difficult for a multi 
parameter fit to be deterministic. However, the above 
demonstrates that our theory offers a sensible fit to the 
data for multiple line cuts in momentum space and the 
temperature dependence of the neutron scattering struc- 
ture factor. In particular, these finite temperature fits 
demonstrate that T* in the range of 40K < T* < lOOif 
and additional thermal fluctuation at 40K compared to 
5K are consistent with the data.. 
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FIG. 5: Single-parameter fitting to momentum dependence of 
INS intensity along a (blue dots) and b (red triangles) axes at 
T = 40K and uj = 3 meV with phonon scatterings. 
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